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Shadow X-ray Magnetic Circular Dichroism Photo-Emission Electron Microscopy (XMCD- 
PEEM) is a recent technique, in which the photon intensity in the shadow of an object lying on a 
surface, may be used to gather information about the three-dimensional magnetization texture in¬ 
side the object. Our purpose here is to lay the basis of a quantitative analysis of this technique. We 
first discuss the principle and implementation of a method to simulate the contrast expected from 
an arbitrary micromagnetic state. Text book examples and successful comparison with experiments 
are then given. Instrumental settings are finally discussed, having an impact on the contrast and 
spatial resolution : photon energy, microscope extraction voltage and plane of focus, microscope 
background level, electric-field related distortion of three-dimensional objects, Fresnel diffraction or 
photon scattering. 


Progress is continuous in the decreasing size and 
increasing complexity of nanosized magnetic systems 
being designed for either fundamental science or de¬ 
vices. Magnetic microscopies are crucial tools to monitor 
and understand the properties of such systems. Vari¬ 
ous types of information are desirable to gather, lead¬ 
ing to multiple criteria to classify microscopies: spa¬ 
tial and time resolution, compatibility with environmen¬ 
tal parameters such as variable temperature and ap¬ 
plied magnetic field, requirements on the sample prepa¬ 
ration and compatibility for ex-situ processing such 
as lithography, correlation with structural information, 
elemental sensitivity, quantity measured (magnetiza¬ 
tion, induction, stray field etc.), sensitivity. The most 
common magnetic microscopies offering spatial resolu¬ 
tion below 50 nm and direct sensitivity to magneti¬ 
zation are X-ray Magnetic Circular Dichroism Photo- 
Emission Electron Microscopy (XMCD-PEEM)[1] and 
(Scanning) Transmission X-ray Microscopy [(S)TXM][2], 
electron holography or Lorentz microscopy [3-5], Scan¬ 
ning Electron Microscopy with Polarization Analy¬ 
sis (SEMPA)[6], Spin-Polarized Low-Energy Electron Mi¬ 
croscopy (SPLEEM)[7-9]. 

Yet another criterion is the volume of the sample 
probed. This criterium is gaining in importance in the 
context of the emergence of three-dimensional (3D) mag¬ 
netic objects and textures. The distribution of magneti¬ 
zation may be truly 3D if along the three directions in 
space the size of a system lies above magnetic characteris¬ 
tic length scales, such as the dipolar exchange length Ad 
for soft magnetic materials, or the anisotropy exchange 
length A u for a hard magnetic material[10]. While this 
is obviously fulfilled in macroscopic materials, the com¬ 
plexity of magnetic textures is such that it cannot be 
measured in detail, and besides it cannot be controlled 
to achieve specific functions. The progress in nanofab¬ 
rication techniques now allows to design suitable sys¬ 


tems, both with top-down and bottom-up approaches. 
Let us give some examples. Flat magnetic elements are 
basic building blocks in spintronics, being patterned with 
great 2D versatility with thin film and lithography tech¬ 
nologies. Large lateral dimensions may give rise to 2D 
magnetic textures, such as the so-called vortex state in 
disks [11]. Such textures are being investigated to de¬ 
sign RF oscillator components, relying on the gyrotropic 
motion of the vortex core driven by a spin-polarized 
current [12]. Stacking several disks provides additional 
degrees of freedom to control the precessional modes or 
the spectrum purity, for instance allowing to consider 
vortices with cores aligned parallel or antiparallel [13]. 
While coupled nanocubes had been imaged by elec¬ 
tron holography in the case of parallel cores [14], stacked 
disks could be imaged recently with magnetic tomog¬ 
raphy holography, revealing fine details in the case of 
(repulsive) antiparallel cores[15]. The topological iden¬ 
tity of such structures with so-called merons has been 
highlighted [16]. Another example is long and cylindri¬ 
cal nanowires and nanotubes [17], which in a dense ar¬ 
ray would be the natural geometry to implement a 3D 
race-track magnetic memory [18]. Simulation and theory 
predicted that a new type of domain wall should exist in 
these systems, with a truly 3D magnetic texture trying 
to close the flux in all three directions[19-21]. Such do¬ 
main walls are expected to display in their core a Bloch 
point, an intriguing magnetic object with a local cancella¬ 
tion of magnetization in the otherwise ferromagnetically- 
ordered material[22, 23]. Accordingly, they are named by 
some[21, 24, 25] Bloch-point wall. These domain walls 
have been predicted to be liable to reach high and ro¬ 
bust velocities [21], due to their specific topology provid¬ 
ing them with a protection against transformations [25]. 
While domains in tubes [26, 27] and conventional trans¬ 
verse walls in wires [28] had been imaged, recently some 
of us provided the experimental proof of the existence of 
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Bloch-point walls using shadow XMCD-PEEM[29] (the 
technique will be introduced in the next paragraph). 
Other examples include magnetization processes inside 
domain walls [30, 31] and dimensional cross-over from vor¬ 
tices to domain walls [32, 33]. 

Let us review again the above-mentioned microscopy 
techniques, in the light of 3D imaging. SPLEEM and 
SEMPA typically probe the topmost atomic layer of mat¬ 
ter. This makes them sensitive to very small amounts of 
material if layered, however hides magnetic information 
in the core of a system. On the reverse, Lorentz, Holog¬ 
raphy and (S)TXM are transmission techniques with a 
penetration depth of the order of 100 nm, providing infor¬ 
mation about volume magnetic textures over this depth. 
However, they have a lower sensitivity as they are not 
strictly surface-based techniques, and average the mea¬ 
sured signal along the path of the beam. Thus, some 
information is lost in the case of magnetic textures vary¬ 
ing along the depth, unless time- and effort-demanding 
tomography procedures are implemented [15]. The prob¬ 
ing depth of XMCD-PEEM and SEMPA is intermedi¬ 
ate, being a few nanometers and related to the mean 
free path of the secondary electrons used for imaging. 
Thus it is not strictly surface sensitive, however not suit¬ 
able a priori to probe magnetic systems in depth. How¬ 
ever, as mentioned above XMCD-PEEM was recently ap¬ 
plied to three-dimensional objects lying on a supporting 
surface[26, 29, 34-36]. As the X-ray beam is tilted with 
respect to the normal to the supporting surface, this pro¬ 
vides magnetic sensitivity both at the surface of the ob¬ 
ject, and gives rise to a shadow on the supporting sur¬ 
face, whose inspection yields information about magne¬ 
tization in the core (Figure la,b). This has been named 
shadow XMCD-PEEM [26]. This provides a technique 
with an interesting hybrid sensitivity, within the set of 
microscopy techniques mentioned above. However, due 
to the three-dimensional shape of the objects considered, 
and the depth- and helicity-dependent absorption of X- 
rays through the structure, the magnetic contrast can¬ 
not simply be interpreted as the projection of magneti¬ 
zation along the direction of the beam, as it is the case for 
the usual surface XMCD-PEEM. For example, Figure la 
shows that the contrast at the surface of a uniformly- 
magnetized 3D object may vary and even change in sign, 
depending on it size. In Ref. [27] the authors simulated 
the contrast in the shadow of a rolled tube. However, 
it was based on an analytical form for the distribution 
of magnetization in a thin sheet, not a simulated micro- 
magnetic configuration. Also, the contrast at the surface 
of the structure was not computed. 

In this manuscript we review specific aspects of shadow 
XMCD-PEEM, and propose a method to analyze the 
resulting images of surface and shadow based on the 
complete micromagnetic structure of an object, to make 
shadow XMCD-PEEM a quantitative technique. The 
manuscript is organized as follows. First the principles 
and implementation of a method to simulate the con¬ 
trast of three-dimensional magnetization textures are de- 


Table I: Absorption coefficients /i of Fe, Ni and Fe 2 oNi 8 o for 
the photon energy set at the L absorption edges of either Fe 
and Ni. Figures for pure elements are derived from Ref.37. 
Figures for the alloy are linear interpolation of figures for 
single elements. 
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scribed. Then we illustrate the simulations with two test 
cases. Comparison with a few experimental cases is then 
made, followed by a discussion of the considerations re¬ 
quired for the quantitative analysis of magnetic contrast 
and spatial resolution. These considerations are largely 
related to finer points of the physics at play, which have 
so far been left aside in the modeling. 


I. METHOD AND IMPLEMENTATION OF THE 
SIMULATIONS 

Our approach consists in considering a given three- 
dimensional magnetization texture in a system, and use 
it as an input to simulate the X-ray absorption spec¬ 
troscopy (XAS, arising from a non-polarized photon 
beam) and XMCD contrasts (the asymmetry ratio aris¬ 
ing from photons with opposite circular polarizations) 
expected at both the surface and in the shadow. The 
magnetic texture may be a simple analytical form, or a re¬ 
alistic distribution of magnetization resulting from a mi- 
cromagnetic simulation. In this section we first describe 
the physical principles considered to convert a magnetic 
texture into a magnetic contrast. Then we detail the 
practical implementation in the numerics. 

A. Principle of the method 

Building a XMCD-PEEM image requires to describe 
mainly three distinct steps including physical and in¬ 
strumental aspects: absorption through matter; photo¬ 
emission of electrons close to surfaces; collection of these 
electrons in the microscope. The way we model each of 
these processes is detailed below. 


1. X-ray absorption 

At any stage when traveling through matter, an X-ray 
beam is associated with a probability of absorption per 







3 





710 


720 730 '' 850 

Photon energy (eV) 


860 


870 


Figure 1: XAS and XMCD with shadows, (a) Shadow 
XMCD (top) on a Fe 2 oNi 8 o cylindrical wire with a modula¬ 
tion of diameter, from 400 nm to 150 nm (bottom: schematic 
of wire, the arrows depicting the direction of magnetization). 
The direction of the beam is indicated by the upper-left wavy 
arrow. Notice the inversion of surface contrast at the back 
side of the wire, when its diameter is the largest (left side) 
(b) Schematic for the dual surface and volume contrast on the 
basis of the test case of magnetization parallel to the X-ray 
beam. The curves below represent the polarization-dependent 
X-ray intensity at the absorption peak as the X-rays propa¬ 
gate through the wire section, (c) Red: absorption spectra on 
the wire across the Fe and Ni L edges, normalized to the back¬ 
ground signal (absorption on the supporting Si surface). Blue: 
inverted and normalized spectra measured in the shadow. 



Figure 2: Collecting the photoelectrons, (a) Illustration 
of electrons escaping the magnetic object. e _ is the secondary 
electron emission direction, and 0 is the angle between the 
imaging axis and the latter direction, (b) Illustration of the 
distortion of the photoelectron trajectory due the electric field 
extracting electrons into the microscope column (see sec.V). 


unit length, /i, determining the mean free path of pho¬ 
tons A = 1 //I. These parameters depend on the compo¬ 
sition of matter, as well as on the photon energy. When 
comparing the model with experiments we will consider 
objects made of permalloy, with composition Fe 2 oNi 80 . 
This will give us the opportunity to discuss the case of 
alloys, and highlight that this gives more freedom for the 
choice of the absorbtion edge. We considered photons at 
the L3 and L2 edges of both Fe and Ni, and used the 
experimentally-determined parameters from the litera¬ 
ture for pure elements and both helicities of photons [37]. 
Based on the very similar volume density of these ele¬ 
ments and alloys, we assumed that the absorption coef¬ 
ficient of permalloy is fi = 0.2/rpe + 0.8/iNi- Absorption 
coefficients at the different edges are summarized in Ta¬ 
ble I. At the Fe L edges, absorption due to Fe is large, 
and the pre-edge absorption of Ni is weak, so that despite 
the low concentration of Fe it is by far the dominating 
contribution. At the Ni L edges the post-edge absorption 
of Fe is in principle no more negligible especially because 
absorption on Ni with a nearly filled 3d band yields a 
moderate intensity. The contribution of Ni is however 
still larger to that of Fe, due to its much larger concen¬ 
tration. 

The progressive absorption of the beam in matter is 
described by integrating its position-sensitive rate of ab¬ 
sorption through each elementary segment with length 
dt 


dlx,a ± 

d£ 


l/i+(l±k-m) + I^_(l=Fk-m) 


Ix,( 


( 1 ) 

where /i + and (i- stand for the absorption coefficients 
for left and right circularly polarized X-rays, respectively. 
This formula takes into account the energy and helicity 
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dependence, in relation with the direction of magneti¬ 
zation in the sample, with k the unit vector along the 
propagation direction. 


2. Local emission of electrons 

We need to estimate the local rate of emission of pho¬ 
toelectrons I e ,a±( r s) at any location r s at the surface, 
resulting from the transmitted X-ray intensity Ix,a±( r s) 
reaching that location, as calculated previously. The ma¬ 
jority of emitted current consists of secondary electrons, 
whose escape depth is only a few nanometers. As this 
length is much smaller than the size of objects of interest 
in shadow-PEEM, and also smaller than any magnetic 
length scale, we used the simplifying assumption that on 
the object / e?cr (r s ) reflects /x,cr( r s) and magnetization at 
the surface, through again the dichroism ratio depending 
on the photon helicity. To the contrary, when impinging 
on the non-magnetic surface, for instance in the shadow, 
the photons give rise to a rate of electrons directly pro¬ 
portional to Tv,cr( r s)- Let us finally discuss the escape of 
electrons from matter. Initially, photo-emitted electrons 
are emitted isotropically and not perpendicular to the lo¬ 
cal sample surface, implying some lateral broadening of 
the electron emission. Thus, we expect that a measured 
image results from the convolution of the signal described 
aforehand, with a function describing these processes. In 
practice however, as the escape depth of electrons is only 
a few nanometers, the expected broadening should not 
exceed these few nanometers, which is much smaller than 
the instrumental resolution (circa 30 nm). Thus, these ef¬ 
fects may be safely neglected. So, at this stage we have 
an estimate / e?cr (r s ) of the local electron emission at each 
point of the nanostructure and its supporting surface. 


3. Intensity on the detector 

We now need to convert the local emission rate / e?cr (r s ), 
into the intensity per unit surface / s (x, ?/), on the detec¬ 
tor (the subscript s standing for screen ). One impor¬ 
tant parameter to make the link between J e?cr (r s ) and 
I s (x,y) is the angular acceptance of the microscope. A 
key parameter is the contrast aperture, whose aim is to 
select electrons escaping the sample essentially along the 
column axis, in order to minimize aberrations. While 
this has no impact for flat surfaces (as for the contrast 
in the shadow), it affects the contrast arising from 3D 
objects. Indeed, the emission of secondary electrons is 
maximum along the normal to the local surface, and the 
total collecting efficiency is lower for tilted surfaces [38]. 
However, the exact angular dependance is not well char¬ 
acterized and may depend sensitively on extraction and 
electron energy, aperture, surface roughness etc. Thus 
we did not attempt to consider a realistic transfer func¬ 
tion, as this would be largely arbitrary. We simply con¬ 
sidered that I s (x,y) = I e , a ( r s ). To understand which 


transfer function this stands for, let us call 0 the angle 
between the microscope axis and the normal to the emit¬ 
ting surface (Figure 2a). The ratio of local sample surface 
over corresponding detector surface scales like 1/cos#, so 
that if all emitted electrons were collected irrespective of 
their direction when escaping the surface, a similar ratio 
would be calculated on the detector. Thus, the choice 
I s (x,y) = / e?cr (r s ) is equivalent to using a cos# collec¬ 
tion function. Notice that while the collection function 
affects the XAS, theoretically XMCD images should not 
depend on this function as they are computed as differ¬ 
ences normalized by the sum. In practice, due to the re¬ 
duced number of photons in areas where the real transfer 
function does not allow the collection of electrons, com¬ 
bined with a background electron level to be discussed 
in sec. IV C, in experimental images the XMCD may be 
sharply decreased in such areas. 

Let us note that the above procedure is not a bijec- 
tion but a surjection, because of the integration along a 
path and also of the projection of magnetization along 
the beam. Thus, one XPEEM image may in principle 
correspond to different magnetic configurations. This, 
along with other issues contributing to image formation 
such as photon scattering, field distortion due to the ob¬ 
ject topography, and background electron intensity, will 
be discussed further below (sec. IV). 

B. Numerical implementation 

Besides analytic test cases, micromagnetic configura¬ 
tions resulting from simulations are used as input to com¬ 
pute the XMCD contrast. For the micromagnetic config¬ 
urations, we use the home-made code FeeLLGood[39]. 
FeeLLGood is based on the finite element method. We 
used material parameters suitable for permalloy: A = 
lOpJ/m, /ioAf s = 1 T. The damping parameter a was 
set to 1 to facilitate convergence, with no impact on 
the results as we only consider states at rest. No mag¬ 
net ocristalline anisotropy was considered. A flow chart 
of experimental and simulation steps is provided in Fig¬ 
ure 3. 

We illustrate here the simulation scheme with cylindri¬ 
cal nanowires, although it can be applied to an arbitrary 
shape. The principle of the numerical method is to con¬ 
sider a ray, and to compute the absorption along this ray 
taking into account the distribution of magnetization in 
the sample. 


1. X-ray absorption 

We first simulate the photon flux Ix, 0 r( r ) with a given 
helicity a at any position in space, inside and outside 
the magnetic system. For this we consider a regular 
planar grid perpendicular to the direction of the pho¬ 
tons of an incident plane wave. From each node a ray 
is launched and intercepts the surface of the cylinder at 
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Figure 3: Flow chart of the analysis method. Key acquisition and simulation steps for the quantitative analysis of shadow 
XMCD-PEEM images. 


two points Pin and P ou t respectively, called the enter¬ 
ing and exit points. In order to calculate the absorption 
along the ray the latter is discretized into segments of 
a given length between P[ n and P out . At any of these 
points we estimate magnetization using an interpolation 
method (see Figure 4). At Pi n and P ou t the magneti¬ 
zation M is interpolated from its values known at the 
nodes of the triangle (the face of the first/last tetra¬ 
hedron) to which it belongs, using a method based on 
the areal coordinates method [40]: M = with 

i = I, II, HI and cq(P in ) = P(P in ,II,III)/P(I,II,III) (Fig¬ 


ure 4c). 5(Pi n , II, III) is the surface of a sub-triangle 

and 5(1,11, III) is the total surface of the triangle, an 
and am are obtained upon circular permutation. For 
points inside the sample, we first determine the cor¬ 
responding element (tetrahedron in the volume of the 
wire, see Figure 4a). The magnetization is then inter¬ 
polated with the method previously described but with 
four nodes and making use of sub-tetrahedrons instead 
of sub-triangles. Then, the progressive absorption along 
the entire path through the magnetic structure can be 
computed. It is done by integrating Eq. (1): Ix,a± 
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Figure 4: Description of interpolation methods. The 

system considered as an example is a cylindrical nanowire, 
(a) presents the method to model the X-ray beam. The red 
line stands for the X-ray beam that crosses the tetrahedron 
elements (in white, from the finite element discretization), (b) 
shows the photon beam, and the wire on the regularly-gridded 
supporting surface, (c) scheme of a triangle to illustrate in 
2D the method of the areal coordinates for interpolating mag¬ 
netization along each ray. 


expj— fp OUt d£ |/i + (l =b k • m) + |/i_(l =F k • m) j. 


2. Local emission of electrons 


The next step is to get the intensity of the photoelec¬ 
trons emitted after the absorption of X-rays. At the sur¬ 
face of the wire, the intensity of the photoelectrons is then 
calculated as the intensity of the photons multiplied by 
a pre-factor and normalized by the density of X-ray flux 
n • k. The pre-factor takes into account the scalar prod¬ 
uct of the local magnetization and the wave vector and 
the local absorption coefficient such that : 


1 


n • k /i+ + /i- \_2 


oM+(l±k-m) 


+ 2^-( 1 T k • m) 


Ix.o 


( 2 ) 


Note that the right hand side is the derivative of X-ray 
intensity along the propagation path [as in Eq. (1)]. The 
physical meaning is clear, because the derivative of Ix,cr is 
indicative of photons absorbed within that incremental 
distance, and photoemission is directly proportional to 
absorption. 


3. Intensity on the detector 

The last step is to infer the electron intensity on the de¬ 
tector on a square grid, from the photoelectron intensity 
previously computed at points of the sample and sup¬ 
porting surface. This is achieved with a linear interpola¬ 
tion, in a similar fashion as described above for P in and 
Pout- Finally, the XMCD-PEEM contrast is computed as 
(J e?cr _ — Ie,a + ) and normalized to the sum (/ e?cr _ +^e,a + )- 
The implementation has been done by using the geome¬ 
try library CGAL[41] for the use of rays, and the nearest 
neighbor searching library ANN [42]. 


II. ILLUSTRATION ON TEST CASES 

In this section we apply the simulation method to two 
test cases of analytical distributions of magnetization: 
transverse uniform magnetization and orthoradial curl¬ 
ing (Figure 5). Curling is a long-existing name is mag¬ 
netism, used to describe an area where the curl of mag¬ 
netization is non-zero, such as around a magnetic vortex, 
or in an object subject to the so-called curling nucleation 
mode as initially introduced [43]. Although they would 
not occur in wires as such, these distributions are chosen 
to illustrate the method, and understand special features 
of contrast which can arise in shadow XMCD-PEEM. 
They are also relevant for the experimentally observed 
magnetic structures observed in cylindrical wires as we 
will show in the following sections. These two situations 
have been described analytically for each point, and we 
checked that an excellent agreement was found with the 
numerical grid method. In Figure 5, the simulated wire 
was suspended above the substrate surface so that the 
complete shadow is collected. Although this may hap¬ 
pen experimentally in some cases (Figure 6b), in most 
cases the wire is in contact with the supporting surface 
so that part of the shadow is not visible on the screen. 


A. Transverse uniform magnetization 

Figure 5a shows ix,o-±( r ) and the resulting dichroic 
absorption for photons going through a wire uniformly 
magnetized along its diameter. The two curves illus¬ 
trate the fact that photons with one of the two polar¬ 
izations is more absorbed than the other, due to mag¬ 
netic dichroism (Figure lb). The dichroic contrast in the 
shadow is therefore opposite to that at the front surface 
of the wire, as it simply reflects the effect in the trans¬ 
mitted photons. Besides, even though photo-emission is 
higher (resp. lower) per photon, for the photons with 
lower (resp. higher) transmission, the dichroism mea¬ 
sured at the back side of the wire may be reversed com¬ 
pared to the front side, for thick-enough wires and thus 
large imbalance after transmission. The dichroic contrast 
at the surface of the wire is illustrated in Figure 5a, con¬ 
stant on the front side and gradually decreasing on the 
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back side. The critical diameter above which the contrast 
reverses may be calculated: 


do = ln(M+/M ~ } . = (3) 

/x + — /i_ \/l — sin 2 0 

where (j) is the incidence angle of the photons (Figure la). 
Obviously, d c depends on the X-ray energy via the ab¬ 
sorption coefficient (i±. In the case of a wire made of 
permalloy, and for a grazing angle (j) = 16° as will be 
reported for experimental comparison, d c is respectively 
70, 140, 20 and 50 nm at the Fe-L3, Fe-L2, Ni-L3 and 
Ni-L2 edges. This explains the dominating inverted con¬ 
trast on the wide-diameter side of the wire in Figure la. 
Note also that the contrast is expected to be larger at 
the center of the shadow than at its border, because the 
length of material probed is larger, and so does the im¬ 
balance of outgoing photons. These facts highlight that 
the contrast does not reflect directly the local direction 
of magnetization, and stresses the need for simulation. 
Practical examples will be provided in sec. IV, dedicated 
to the analysis of contrasts. 


B. Orthoradial curling 

The case of orthoradial curling of magnetization (see 
Figure 5b) is directly relevant for one type of domain wall 
in cylindrical wires: the Bloch point wall [25, 29]. At the 
bottom part of the wire, magnetization is mostly pointing 
left, while at the top part it is mostly pointing right. This 
leads to opposite contrasts on either side of the shadow. 
The center of the shadow has no XMCD contrast, as at 
all points through the wire diameter the X-ray direction 
is perpendicular to the magnetization direction. This 
dipolar contrast is a clear signature of orthoradial curl¬ 
ing. To the contrary, the dichroic contrast at the surface 
of the wire is maximum close to its top, where the beam 
is tangent to its surface. It decays on both sides, with a 
slight negative value on the front side due to the tilted 
incidence of the photons, and a possible inversion of con¬ 
trast on the backside depending on the total absorption. 
Thus, the contrast is largely monopolar as in the case 
of uniform transverse magnetization, which could for in¬ 
stance naively be expected from a transverse wall with 
the transverse component aligned with the beam direc¬ 
tion. Thus, ascribing the surface contrast to a trans¬ 
verse wall or a Bloch-point wall may remain ambiguous. 
This example shows that inspection of the shadow may 
be crucial to get information about a three-dimensional 
configuration of magnetization. 


III. COMPARISON WITH EXPERIMENTS 

To put the model into practice, we take the example of 
experiments using shadow XMCD-PEEM on cylindrical 
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Figure 5: Illustration of XMCD-PEEM post¬ 

processing on two test cases. The first test case is (a) 
uniform magnetization across the wire and parallel to the X- 
ray beam; the second case (b) is orthoradial curling. The 
lower part presents the photon density and XMCD for each 
case at the surface of the wire (yellow background) and in the 
shadow (pink background). The modeled wire is suspended 
above the surface, so that the entire shadow is visible. Note 
that the lateral scale is expanded by a factor sin(16°) ^3.6 
in the shadow, thanks to the grazing incidence. 


wires, consistent with the test cases discussed above. In 
this section, we describe the experimental setup, and dis¬ 
cuss the correspondence of the experimentally observed 
images with those from the simulations. 


A. Experimental details 

The samples considered are permalloy cylindrical wires 
electroplated in self-organized anodized alumina tem¬ 
plates. The alumina matrix is dissolved and wires are 
dispersed on a naturally-oxidized Si supporting surface. 
The wires are aligned along a preferential in-plane direc¬ 
tion thanks to an in-plane magnetic field applied during 
dispersion. Their diameter, possibly modulated along 
the length, ranges from 50 nm to several hundreds of nm. 
The length of the wires is typically a few micrometers [29]. 
Element-selective XAS and XMCD-PEEM were carried 
out at the spectroscopic photoemission and low-energy 
electron microscope [44] operated at the undulator beam¬ 
line Nanospectroscopy of Elettra, Sincrotrone Trieste. 
The photons impinge on the surface with a grazing angle 
(p = 16 °. Spectroscopy was performed across the L edges 
of either Ni or Fe, using elliptically-polarized radiation as 


































a probe. Series of several tens of images with an exposure 
time of few seconds are recorded, drift-corrected and fi¬ 
nally co-added. This yields a high signal-over-noise ratio 
while limiting drift effects, providing images with a spa¬ 
tial resolution of the order of 30 nm. The level of circular 
polarization at the Fe and Ni L edges was estimated to be 
around 75 % as the X-ray beam is produced by a higher 
harmonic of the undulator source. 


B. Experimental test cases: curling structures 

Two types of domain walls may be expected in cylin¬ 
drical nanowires: of mixed transverse-vortex type for 
diameter below typically 7A^ (Ad = y/2A//ioM s 2 ), 

and of Bloch-point type for diameter above typically 
7Ad [21, 25]. The former is reminiscent of domain walls 
already known in flat strips [45, 46], while the latter is 
specific to wires with large dimensions. Consistently, 
we observed two types of contrast for domain walls in 
nanowires, which we could ascribe to these walls [29]. 
Here we illustrate the shadow technique with the Bloch- 
point wall (Figure 6). 

Figure 6a shows the XMCD-PEEM image from a 
magnetic wire of diameter 90 nm, suspended above the 
substrate in the imaged region. A good agreement is 
found with the contrast delivered by our model (Fig¬ 
ure 6b), based on the simulated micromagnetic config¬ 
uration (Figure 6c). Thus, this domain wall can be un¬ 
ambiguously ascribed to a Bloch-point wall. The strik¬ 
ing feature is the dual bright and dark contrasts in the 
shadow, revealing a local orthoradial curling as already 
seen in Figure 5b. The symmetry with respect to a 
plane perpendicular to the wire axis shows that curling 
is purely orthoradial, which is compatible only with the 
Bloch-point wall [25]. A quantitative comparison may be 
made with a cross-section (Figure 6d). The contrast has 
been normalized so that the maxima coincide at the sur¬ 
face of the wire. Compared to simulations, the experi¬ 
mental cross-section is wider by approximately 25 nm on 
either side, which is consistent with the expected instru¬ 
mental spatial broadening. The agreement is however 
excellent at the surface of the wire, especially the rather 
sharp maximum and its location away from the central 
part and towards the backside of the wire (Figure 6d, 
insert). This feature is explained as follows. The XAS 
should have the shape shown on the top part of Figure 5b. 
Upon convolution with the resolution function, this ini¬ 
tially asymmetric XMCD shape gives rise to a maximum 
displaced towards the backside of the wire. Note also 
that in the shadow, the cross-section is clearly antisym¬ 
metric, as expected. The cancellation of contrast at the 
core of the shadow should coincide with the location of 
the Bloch point. The experimental contrast is however 
lower than expected in theory, which will be discussed in 
sec. IV. Of importance is the fact that the lateral scale 
of the wire is expanded by a factor 1/sin(16°) ~ 3.6 in 
the shadow, thanks to the grazing incidence. In princi- 
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Figure 6: Comparison of experiment and simulation for 
the Bloch-Point wall, (a) Experimental and (b) simulated 
contrasts at the surface and in the shadow (resp. left and right 
in both images) at the Fe L3 edge for a wire of diameter 90 nm. 
This diameter has been determined from the XAS of the wire 
shadow, (c) Micromagnetic simulation of a Bloch-point wall, 
used as an input for simulating the contrast (d) Cross-sections 
for (a) and (b). The insert shows the XAS (dotted line) and 
XMCD (full line) spectra for the experimental contrasts (a). 


pie this promises an increase of spatial resolution of 3.6 
along one direction, however issues of signal-over-noise 
ratio may limit this gain, which will be addressed in the 
next section. 

Besides domain walls, three-dimensional non-uniform 
distributions of magnetization are also expected at diam¬ 
eter modulations and at the ends of wires, driven by the 
reduction of magnetostatic energy. Curling of magneti¬ 
zation around the wire axis has been predicted at such 
locations [47, 48], however remaining elusive experimen¬ 
tally so far. What has been reported are hints for the 
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Table II: The absorption coefficient values used for the simu¬ 
lations in Figure 7, (k-n). 



(k,l) 

(m,n) 

/i_ (nm 1 ) 

0.018 

0.036 

M+ (nm -1 ) 

0.01 

0.02 


spread of charges, which however could also be argued 
to take the form of other flux-closing structures [49, 50]. 
Shadow PEEM again provides a direct proof for the exis¬ 
tence of such buried structures. Figure 7(c-f) shows the 
XMCD contrast of a wire at various absorption edges. A 
close-up view of the end of the wire is displayed in Fig¬ 
ure 7b, along with a simulation derived from a curling 
end domain. Thanks to the comparison we can formally 
identify the contrast at the end as arising from an ortho- 
radial curling structure. 


IV. DISCUSSION ON CONTRAST 

In this section we discuss in more detail several instru¬ 
mental aspects specific to the shadow imaging geometry, 
which have an impact on the magnetic contrast or spa¬ 
tial resolution. Of special importance for the shadow 
imaging of 3D objects are the plane of focus, the start 
voltage (STV, a voltage bias applied to the sample and 
which determines the electron kinetic energy), the mi¬ 
croscope background level, the photon energy, and the 
Fresnel diffraction of X-rays. 


A. Microscope settings 

While a rather flat surface may be entirely set close 
in focus, the case of three-dimensional objects lying on a 
surface is different. The depth of focus of the instrument 
is several micrometers, large enough so that the top of 
the wire and the supporting surface may both be in focus. 
In practice however this could not be achieved, which we 
understand as resulting from the wire curvature acting 
as a lens. For each image, one may thus decide to set the 
focus anywhere between the top surface of the wire, and 
the supporting surface. For instance, setting the focus on 
the supporting surface has a dramatic effect on loosing 
sharpness and therefore XAS and XMCD contrast on the 
wire, due to its small lateral size (Figure 8a, b). Blurring 
effects are decreased upon increasing the start voltage, 
for reasons described below. 

Second, electrons are extracted into the imaging col¬ 
umn with a voltage 18keV— STV. The start voltage 
(STV) is an additional bias, which is related to the elec¬ 
tron kinetic energy (with an additional offset due to work 
function difference between sample and the LaBg source 
used as reference for the energy scale). The non-planar 
wires we use create a non-uniform potential profile of the 


a 1 > /A " A A ^ 








500nm 





Figure 7: Role of the absorption coefficient on the con¬ 
trast. (a) LEEM image of a wire with diameter 120 nm (b) 
left end of the wire, with both experimental and simulated 
XMCD contrast at the Fe L3 edge, the latter revealing a curl¬ 
ing micromagnetic configuration, (c-f) experimental XMCD 
contrast of this wire, (g-j) are the XMCD contrast com¬ 
puted from the same experimental XAS, however from which 
the background level has been removed. Simulated XPEEM 
images (k,m) and their profiles (l,n) for a 90 nm-diameter 
wire with a Bloch point wall and different absorbtion coef¬ 
ficients (see Table II). For (c,d,g,h) the contrast is 6%, it is 
9% for (f,j) and 5% for (e,i). 
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surface, which distorts the emitted electron wave. The 
lower the electron energy the more pronounced this dis¬ 
tortion is. That may explain why the image quality both 
in the wire and in the shadow is better in Figure 8e than 
Figure 8d. 

A fine tuning of the start voltage may also be used to 
enhance the signal originating from either the wire sur¬ 
face or the shadow. Indeed the materials giving rise to 
photoemission are different (here permalloy for the mag¬ 
netic object, and Si for the supporting surface), as well 
as their capping, so the energy distribution and yield of 
secondary photoelectrons are different (Figure 8f). In 
the present case a lower start voltage (~ 2.0 eV) maxi¬ 
mized the number of electrons emitted from the shadow, 
whereas a higher start voltage (^ 2.8 eV) resulted in 
a higher intensity emitted from the wire surface (Fig¬ 
ure 8e). XAS being the measure of emitted electrons, 
these effects of start voltage are directly transferred to 
the XAS image. As XMCD is a difference normalized 
to a sum, its magnitude should not depend on the num¬ 
ber of electrons and thus be insensitive to the choice of 
start voltage. However, as will be argued in sec. IV C, in 
practice a lower number of emitted electrons reduces the 
XMCD signal with respect to the computed value, so the 
start voltage also has an impact on the relative level of 
contrast on the wire versus the shadow. 




B. Photon energy 

For 3d ferromagnetic metals the photon energy needs 
to be tuned close to the maximum of an L3 edge to maxi¬ 
mize dichroism. Attention should be paid to the fact that 
event a slight surface oxidation, induces a sub-structure 
in the L3 peak. This slightly shifts the XAS maximum 
with respect to a metallic spectrum (Figure lc), while 
XMCD remains maximum at the location of the metal 
peak. Accordingly, we found more precise to use an ab¬ 
sorption spectrum taken in the shadow to set the working 
photon energy, as this probes the bulk of the wire, with 
no oxidation. Notice that this choice also maximizes the 
dichroic contrast at the surface of slightly-oxidized wires. 
In practice, we worked at the Fe L3 edge, which despite 
the low concentration of Fe, proved to yield a larger con¬ 
trast than the Ni L3 edge, both at the surface and in the 
shadow of the wires. This will be addressed in sec. IV C. 

We showed in sec. II that a positive XMCD contrast 
on the wire should be associated with a negative XMCD 
contrast in the shadow: a larger absorption and thus 
loss of photons of a given helicity is associated with an 
enhanced number of emitted electrons. The number of 
photons going through the magnetic object depends on 
the dimensionless quantity d/i±, with d its the depth, 
here the diameter of the wire. The effect of varying this 
quantity is illustrated by a movie of the wire and its 
shadow upon ramping the photon energy from below Fe 
L edges to above Ni L edges (see supplementary mate¬ 
rials). For large enough d/i±, the imbalance of photons 


Figure 8: Role of the plane of focus and the start volt¬ 
age. XMCD contrast taken at the Fe L3 edge and with start 
voltage 2.1 eV, with focus (a) on the surface and (b) on the 
wire, (c) is the XAS, superimposed with wire and shadow 
areas as a guide to the eye. (d) XMCD with a focus on the 
surface, and start voltages 2.0 eV and (e) 2.8 eV. (f) XAS on 
the wire and the shadow versus the start voltage. 
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ratio (see discussion in sec.VB). 


C. Background level in PEEM imaging 

On a theoretical basis the contrast in the shadow could 
reach arbitrarily high values for high /id, which however 
comes at the expense of much reduced intensity. This 
is the principle of some polarizers, for example for the 
helicity of X-rays [51] or spin of electrons [52]. Aside 
from obvious issues arising from the signal-over-noise ra¬ 
tio, we found out that in practice an instrumental effect 
limits the contrast. Figure 9 shows the level of XAS 
on a broad wire and its shadow. For this broad wire 
fi±d ^ 1, so that the intensity in the shadow should 
be vanishingly small. To the contrary, although the in¬ 
tensity reaches a plateau inside the shadow, it remains 
close to 7 % of the intensity over the free supporting sur¬ 
face (Figure 9d). This intensity is not related to the 
background electronic level of the camera, which is al¬ 
ready subtracted from the images. Instead, it reflects 
electrons that truly impinge on the detector. The phys¬ 
ical origin of this background is not straightforward, as 
it was found to be only weakly affected by changing set¬ 
tings of the LEEM. In particular the field-of-view aper¬ 
ture was found to be unrelated, although rejecting elec¬ 
trons from the imaging column arising from outside the 
field of view, to avoid their incoherent contribution to 
the image. The contrast aperture, affecting the angular 
collection of the microscope, did not have a sizeable im¬ 
pact either. Thus although its origin is not clear, an (a 
priori) helicity-independent background intensity reduces 
the computed XMCD as its difference in the numerator 
is zero, while its sum in denominator is non-zero. If the 
background intensity / e? b is known, then a more accu¬ 
rate view of the true XMCD is achieved by computing: 
^XMCD,0 = (^e,cr_ -^e,a + )/(^e,cr_+^ e ,a + -2/ e ,b)- This has 
been done in Figure 7(g-j). It is striking that the con¬ 
trast in the shadow is enhanced, as expected from theory. 
These comparisons also illustrate that working at the Fe 
edges yields in practice a higher contrast than at the Ni 
edges, whereas a similar contrast would be expected for 
permalloy as computed from the tabulated absorption 
coefficients (Table I). This is explained simply by the ex¬ 
istence of the background level. 


D. Scattering effects 

In general, interaction of X-rays with matter can be de¬ 
scribed via the complex atomic scattering factors. The 
real and imaginary parts give rise to the Faraday rota¬ 
tion of the photon polarization vector and to magnetic 
dichroism, respectively, as the photon beam propagates 
through the magnetic material. The two are related by 
the Kramers-Kronig transformation, and they are com¬ 
parable in magnitude at the Fe L3 edge [53]. 





Path length (nm) 


Figure 9: Scattering effects. Wire image of the intensity 
for (a) the XMCD and (b) the XAS. (c) shows the wire and 
shadow zones on the XAS images. The plot (d) presents the 
intensity at the cross-sections shown in (a) and (b). Dotted 
lines are guides to the eyes and limit wire (left) and shadow 
areas (right). 


Until now, we have considered only the X-ray absorp¬ 
tion coefficient, which is proportional to the imaginary 
part of the forward scattering amplitude via the optical 
theorem [54]. Instead, as we noted above, a proper treat¬ 
ment should include the full scattering process. Indeed, 
intensity oscillations near the shadow edge are visible in 
Figure 9 due to Fresnel diffraction from the wire. Fur¬ 
thermore, the Fresnel fringes also show a dichroic signal. 
The sign of this dichroic signal is opposite to that ob¬ 
served within the shadow, as expected from the inverted 
absorption signal in transmission. 

Nevertheless, the shadow (or the substrate) is in the 
very near field of the wire, and coherent scattering effects 
are limited to the shadow edge. The Fraunhofer region 
at this wavelength (about 1.8 nm) and for a wire diame¬ 
ter below 100 nm does not before a few tens of microns. 
Therefore, our analysis relating the dichroism within the 
shadow to the absorption coefficient is valid except at the 
very edge of the shadow [55]. 


V. DISCUSSION ON SPATIAL RESOLUTION 

In this last section we discuss spatial resolution effects 
specific to shadow imaging. 
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Figure 10: Signal-over-noise ratio. Schematic to describe 
the notations used in the calculation (see text). 


A. Electric-field related distortion 

First, it is known that LEEM images of non-planar sur¬ 
face are distorted [1]. The physical phenomenon is that 
secondary electrons escape the material perpendicular to 
the local surface, on the average. Thus their trajectory 
is curved through the extraction electric field, the curva¬ 
ture of the object providing the same effect as a lens (Fig¬ 
ure 2b). Besides curvature, the complex topography of 
the sample acts as a cathode for the accelerating voltage. 
Therefore the accelerating voltage is not uniform across 
the surface. This creates a significant distortion to the 
outgoing low-energy electron wave. The LEEM image 
of convex and topographically-complex objects such as 
the wires considered here, is therefore expected to dis¬ 
play a larger size than the real object. Obviously this 
phenomenon is largely absent in the shadow, as the tra¬ 
jectory of photons is only weakly affected by the circular 
shape of the wires, due to the optical index being very 
close to unity. For wires lying perfectly on the supporting 
surface this effect cannot be checked, because the direct 
and shadow contrasts overlap. In a few cases a gap was 
found between a wire and the surface, large enough to 
separate the direct and shadow areas (Figure 6b). For 
such cases the apparent width of the wire deduced from 
XAS images was indeed about 50% larger on the wire 
than in the shadow. To minimize systematic errors, the 
figures for wire diameter mentioned in the manuscript 
and used for simulations, have always been those deduced 
from the shadow and deconvolved from the expected 
30 nm experimental resolution. 


B. Signal-over-noise ratio in the shadow 

Not only is the shadow more faithful as just discussed, 
but it may promise for an increase of spatial resolution by 
a factor 1/sin(16°) ~ 3.6 thanks to the projection with 
a rather grazing incidence. This would bring the spatial 
resolution along one direction below 10 nm. A practical 
limitation for this gain is the lower number of electrons 
collected in the shadow, degrading the signal-over-noise 
ratio as estimated below. Let N be the number of elec¬ 
trons emitted from the supporting surface under direct il¬ 


lumination, per given time and area. TV b and N s h are sim¬ 
ilarly the number of electrons contributing to the back¬ 
ground level, and those contributing to the shadow and 
related to photons transmitted through the wire, again 
per unit area and time. In our case 7V b /7V « 0.07 and 
N sh /N ~ exp (—/id), with d the diameter of the wire. 
The shortest spatial variation that can be expected on 
the detector is the instrumental resolution cri nstr , result¬ 
ing in a slope f = tA s h/cri nstr with t the averaging time. 
When analyzing experimental data, the possible error on 
lateral resolution a x resulting from the vertical error bar 
(j y is such that Gy/a x = f' (see Figure 10). We thus 
have: a x = cr instr (cr y /£A s h). Taking into account that 
a y = yftNb + \/L/Vsh, one finally gets: 

'■-^( i+ TS) ,4) 

In the absence of background level, Eq. (4) boils down 
to the usual statistics: a x = cr in str/VL/V s h. Thus in 
theory an image of quality similar to that outside the 
shadow with integration time to, could be obtained at 
the expense of an increase in integration time up to t s h 
such that \/4hTV s h = y/toN, so with an increase of time 
N/N s h = exp (/id). This ratio is of the order of 10 2 — 10 4 
for a wire with a diameter of 100 nm at the Fe-L3 edge. 
However in the case of non-zero background level, fluc¬ 
tuations in 7V b also contribute to the increase of a y and 
thus of a x . If A s h becomes small compared to 7V b , based 
on Eq. (4) the time required to get an image of similar 
quality is £ s h = £o(A/A s h) 2 (A b /TV), thus with now the 
power law exp(2/id). For a diameter 100nm, the power 
law is now proportional to 10 4 — 10 8 at the Fe L3 edge 
which becomes prohibitively large. It is the same effect 
of limited statistics, which limits the signal-over-noise ra¬ 
tio of dichroism in the shadow. While in sec. IV B we saw 
that dichroism should asymptotically reach 100% in deep 
shadows (/id >> 1), one can show that in practice the 
signal-over-noise is maximum for /id ~ 1. 

VI. DISCUSSION AND CONCLUSION 

At this stage a comparison may be made with the 
transmission X-ray microscope (TXM). Indeed both 
TXM and shadow-PEEM allow for probing the volume 
magnetization integrated along the photon beam. An ad¬ 
vantage of TXM is its all-photon basis, making it easily 
compatible with applied magnetic fields. Also, the sam¬ 
ple may be rotated to some extent, gaining information 
on different directions of magnetization or integration. 
On the reverse, shadow-PEEM provides the combination 
of surface and volume information, which may be crucial 
to solve complex three dimensional magnetization distri¬ 
butions. The potential increase of spatial resolution is 
also unique. Experiments may even be designed with 
magnetic objects tilted on purpose to a chosen angle to 
make the best use of this gain. 
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To conclude, we have discussed quantitatively physi¬ 
cal and instrumental features specific to shadow-PEEM 
imaging of three-dimensional objects lying on a support¬ 
ing surface. We have considered in more detail XMCD 
imaging and simulation of the expected contrast from 
micromagnetic simulations. This technique uniquely pro¬ 
vides the combination of surface and volume sensitivity in 
the signal measurement, with an enhanced XMCD con¬ 
trast and several-fold gain in spatial resolution along the 
beam direction for the latter. Several effects mentioned 
need however to be considered to extract true spatial and 
contrast information such as plane of focus, extraction 
voltage, electric field distortion and electron background 
level. While we illustrated the method with experiments 
we performed on cylindrical nanowires, it can be ap¬ 
plied to any object of size 10—200 nm imaged by shadow 
XMCD-PEEM, such as those already reported[26, 34- 
36] , provided that material parameters such as exchange 


stiffness and magnetization are known, to perform reli¬ 
able micromagnet ic simulations. In particular, sufficient 
geometrical information about the sample remains re¬ 
quired, as its projected shadow does not characterize fully 
its shape. 
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